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I.  INTRODUCTION 


A.  MOTIVATION 

Binary  computing  systems  have  been  used  for  about  30  years.  The  very-large-scale- 
integration  (VLSI)  technology  has  resulted  in  increasingly  smaller  circuits.  However, 
the  pace  of  reduction  has  showed  because  of  two  major  problems,  bus  connection  and 
pin  limitation.  In  recent  years,  multiple-valued  logic  has  been  proposed  to  solve  those 
problems  [Ref.  1,  2],  Since  the  programmable  logic  array  (PLA)  is  a  basic  tool  for 
binary  VLSI  design,  it  is  expected  to  be  important  in  multiple- valued  logic  (MVL).  This 
observation  has  inspired  the  work  reported  in  this  thesis. 

The  truncated  sum  (TSUM)  operation  of  multiple-valued  logic  is  easily  imple¬ 
mented  in  the  charge-coupled  devices  (CCD)  technology  [Ref.  3.  4].  Several  heuristics 
have  been  developed  for  the  multiple- valued  logic  minimization  involving  TSUM  and  each 
claims  some  advantages  in  specific  examples,  but  none  of  them  is  consistently  better  than 
the  others  [Ref.  5,  6,  7,  8,  9,  10].  Heuristic  methods  are  interesting  because  of  exact  min¬ 
imization  methods  are  extremely  time-consuming.  Until  now,  only  Tirumalai  and  Butler 
[Ref.  9]  have  analyzed  the  performance  of  different  MVL  minimization  algorithms.  Their 
study  is  limited  to  two- variable  four-valued  functions.  Three  open  questions  exist: 

1.  In  the  case  of  two- variable  four-valued  functions,  can  we  design  an  algorithm  which 
is  better  than  previously  proposed  algorithms? 

,2.  What  is  the  algorithm  performance  if  we  increase  the  number  of  variables  (greater 
than  two)? 

3.  What  are  the  computation  time  requirements  for  the  various  algorithms? 

With  the  computer  software  developed  at  Naval  Postgraduate  School  called  HA.M- 
LET  (Heuristic  Analyzer  for  Multiple- valued  Logic  Expression  Translation)  [Ref.  10], 
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users  can  easily  investigate  their  own  heuristics.  For  example,  the  neighborhood  de¬ 
coupling  (ND)  algorithm  described  in  this  thesis  is  built  as  one  independent  option  of 
HAMLET.  This  thesis  reexamines  the  algorithms  used  in  HAMLET  and  analyzes  the 
new  method,  the  ND  algorithm.  It  also  investigates  each  algorithm  under  four  different 
function  settings  and  summarizes  their  performance  and  computation  time. 

B.  BACKGROUND 

The  heuristic  approaches  in  HAMLET  can  all  be  classified  as  greedy  algorithms. 
The  HAMLET  execution  procedure  of  these  algorithms  is  abstracted  as  follows.  Formal 
definitions  will  be  covered  in  the  next  chapter.  Let  /  be  a  multiple-valued  function  and 
let  a,  a  minterm,  be  an  assignment  of  values  to  the  variables  of  /  such  that  /  ^  0. 

j*  XxxXxxXXxXX***^*******X*.*X^*X****X*****ic.ifX***XX*X**XXXX**xxx*xxxxxxxxxxx 

Input:  let  the  M  be  the  set  of  minterms  of  a  function  f; 

Output:  the  minimized  sum  of  product,  S,  of  the  original  function: 

************************************************************************  */ 

{ 

S  +-  <j>. 

While  ( M  ±  <f>)  do  { 

pick  one  minterm  a  from  A/; 
find  an  implicant  Ia  which  covers  a; 

S  <—  7a  U  S; 
subtract  Ia  from  /; 

} 

} 
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TABLE  1.1:  Summary  of  Three  Heuristic  Algorithms 


Heuristic  Algorithm 

Choice  of  Minterm 

Choice  of  Implicant 

Pomper  and  Armstrong  [Ref.  5] 
(1981) 

Random 

Drives  Most  Minterms  to 

0  or  don’t-care 

Besslich  [Ref.  6] 

(1986) 

Smallest  Weight 
(Most  Isolated) 

Drives  Most  Minterms  to 

0  or  don’t-care 

Dueck  and  Miller  [Ref.  7,  8] 
(1988) 

Largest  IF 
(Most  Isolated) 

Largest  BCR 

Table  1.1  shows  three  previously  proposed  algorithms.  Each  algorithm  differs  from 
the  others  in  the  manner  of  picking  the  minterms  (a)  and  finding  the  implicants  (Ia). 
For  example,  in  the  above  recursive  search  procedure,  the  input  function  expression  /  is 
evaluated  at  minterm  a.  Next  an  implicant  Ia  is  chosen  which  covers  a.  Then,  implicant 
Ia  is  added  to  output  solution  set  5,  and  is  subtracted  from  function  /. 

Pomper  and  Armstrong  introduced  a  heuristic  algorithm  that  picks  a  randomly  (as 
long  as  a  is  in  the  set  of  minterm  M )  and  finds  IQ  (as  long  as  Ia  covers  a)  w'hich  drives  the 
most  minterms  to  0  or  don’t-care  when  Ia  subtracted  from  function  /  [Ref.  5].  In  1986, 
Besslich  presented  a  direct  cover  algorithm,  according  to  weight  transformations.  The 
Besslich  algorithm  picks  a  with  the  smallest  weight  (most  isolated  minterm)  and  finds  Ia 
which  has  a  lowest  cost  per  minterm  covered  (i.e.,  which  drives  the  most  minterms  to  0 
or  don’t-care)  [Ref.  6].  In  1988,  Dueck  and  Miller  presented  another  algorithm  that  picks 
a  from  M  if  a  has  the  highest  isolated  factor  (IF)  and  then  finds  the  Ia  which  directly 
covers  a  such  that  the  break  count  reduction  (BCR)  is  maximum  [Ref.  7,  8]. 

Each  algorithm  has  some  advantages  in  specific  examples  [Ref.  9].  The  ND  algo¬ 
rithm  is  characterized  by  adopting  the  advantage  from  each  algorithm  and  fully  utilizing 
the  properties  of  the  truncated  sum.  The  ND  algorithm  is  an  improvement  to  the  Dueck 
and  Miller  algorithm  with  revised  decision  rules  for  making  selections  of  minterms  and 
implicants. 
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C.  THESIS  OUTLINE 

A  summary  of  MVL  definitions  for  truncated  sum  minimization  are  introduced  in 
Chapter  II.  The  notations  and  definitions  of  Chapter  II  also  help  us  in  explaining  the 
algorithms  appearing  in  subsequent  chapters.  The  neighborhood  decoupling  algorithm 
is  presented  in  Chapter  III.  Chapter  IV  and  V  discuss  the  performance  and  computa¬ 
tion  rimes  of  the  neighborhood  decoupling  algorithm  with  Pomper  and  Armstrong’s  and 
Dueck  and  Miller’s  algorithm  for  different  function  parameters. 
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II.  NOTATIONS  AND  DEFINITIONS 

Although  multiple-valued  logic  minimization  for  truncated  sum  has  been  discussed 
in  resent  years  [Ref.  7,  8,  9,  10,  11],  a  distinct  and  complete  formal  definition  for  truncated 
sum  MVL  minimization  is  not  available.  In  view  of  this,  we  studied  the  truncated  sum 
MVL  minimization,  then  summarized  and  give  several  definitions  that  will  be  presented 
in  Section  A.  It  will  also  help  us  to  explain  the  algorithm  given  in  the  next  chapter.  We 
will  report  our  observations  on  truncated  sum  properties  in  Section  B,  and  the  definitions 
for  ND  algorithm  will  be  discussed  in  Section  C. 

A.  DEFINITIONS  FOR  TRUNCATED  SUM 
Definition  1: 

Let  X  =  {  x\,  x2, . . . ,  xn  }  be  a  set  of  n  input  variables  where  x,  takes  on  values 
from  1Z  =  {0, 1, . . . ,  r  —  1}.  An  n- variable  r- valued  function  /  is  a  mapping 

f  -.IT  — >ftU{r}.  [Ref.  9] 

Here,  r  is  a  don't  —  care  value;  it  can  be  chosen  freely  from  any  of  the  logic  values, 

0, 1,  •  •  • ,  r  —  1. 

Definition  2:  MIN 

The  MIN  [Ref.  9]  function  is  denoted  as  f(x  1,12)  =  xix2  which  evaluates  to  the 
minimum  value  of  its  arguments.  For  example,  if  H  =  {0,1, 2, 3},  then  /(1,2)  =  1  and 
/( 0, 3)  =  0.  A  minterm  is  an  assignment  of  values  to  xj,  x2, . . . ,  xn  such  that  /(x)  ^  0. 
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Definition  3:  Literal 


The  literal  operation  of  a  variable  x  is  defined  as: 

v  =  | 

Definition  4:  Truncated  Sum  (TSUM) 


r  —  1  if  a  <  x  <  b 
0  otherwise.  [Ref.  11] 


(2.1) 


The  truncated  sum  (TSUM)  operation  is  defined  as: 

TSUM{ ii,x2)  =  xx  +  x2  =  min(xi  +  x2,r  -  1).  [Re/.  11]  (2.2) 

The  two  +  signs  in  this  expression  are  different.  The  leftmost  denotes  the  TSUM 
operation,  while  the  rightmost  denotes  ordinary  add  two  of  logic  values  which  are  viewed 
as  integer.  The  TSUM  obeys  the  associative  and  commutative  rules.  For  example,  if  IZ 
=  {0,1, 2, 3},  then  TSUM(1,2)=3  and  TSUM(2,2)=3. 

These  definitions  are  inspired  by  the  fact  that  CCD  implementation  supports  TSUM 
naturally  [Ref.  3,  4], 

Example  1: 

For  example,  lx\  is  a  literal  and  takes  value  of  3  when  1  <  xx  <  3.  However, 
function  2  1x?  takes  a  value  of  2  based  on  to  the  definition  of  MIN.  Similarly,  the 
product  term  2  3x l  °X2  is  a  function  that  is  2,  when  xi  =  3  and  x2  =  0. 

Definition  5:  Product  Term 


A  product  term  p  is  the  MIN  of  one  nonzero  constant  c  €  R,  and  one  or  more 
literal  functions.  In  general,  a  product  term  is  defined  as: 

1  <  k  <n.  <2'3> 


p  =  C  HXiJ 


The  constant  or  coefficient  c,  in  a  product  term,  effectively  scales  the  term.  For 
each  variable  x,,  we  say  the  window  size  of  the  literal  is  jk  —  **  +  1.  We  use  the 

terms  product  term  and  implicant  interchangeably  in  this  thesis. 
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Definiton  6:  Sum-of-Products  Expression 

A  sum-of-products  expression  is  p\  +  p2  -f  • . .  +  Pn  for  some  integer  N,  where  p,  is 
a  product  term. 

Definition  7:  Minterm 

A  minterm  a  is  of  the  form  c  aiXj'  02 x 03 . .  .  anx“n  where  a,  €  71  and  constant 
c  €  1Z—  {  0  }.  We  say  the  coordinate  of  q  is  <  ai,a2, . . .  an  >•  We  denote  the  value  of 
minterm  a,  g{a),  as  the  nonzero  constant  c. 

A  product  termp  =  c  ‘‘x^1  ,2x2J2  •  •  ■'nxJn  can  be  decomposed  into  n*=i  (,7fc  —  U  +  l) 
minterms.  We  say  p  generates  those  minterms.  Given  a  product  term  p,  the  set  of 
minterms  generated  from  p  is  denoted  by  MSP.  If  the  number  of  elements  in  MSPl  is 
greater  than  that  in  MSP 2,  we  say  pi  covers  a  larger  area  than  p2-  Given  a  function  /, 
the  set  of  minterms  generated  from  its  product  terms  is  denoted  by  MS/. 

Definition  8:  Saturated  Minterms  (SAT) 

Given  a  minterm  a  generated  from  the  original  function  to  be  minimized,  if  g(a)  = 
r  —  1,  then  a  is  a  saturated  minterm.  Let  SAT  be  the  set  of  all  saturated  minterms  of  a 
function. 


Example  2: 

If  the  input  function  to  be  minimized  is  expressed  as  follows, 

/  =  3  °Zj  xx\  +  2  lx\ °x°  +  3  2xf  +  2  2ii  2x§  +  1  2x2  3x\  +  1  °x°  2x2, 

the  MSf  can  be  represented  as  11  minterms  in  Figure  2.1.  We  mark  a  saturated  minterm 
with  a  dot  in  the  figure. 
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Figure  2.1:  Map  for  Example  2,  3,  4;  Step  1  of  Table  3.2 

Lemma  1  Given  a  minterm  a  the  maximum,  number  of  implicants  which  covers  a  is 

0(r2n). 

Proof:  Consider  a  variable  (axis)  X{  of  a.  Any  implicant  (Ia)  that  covers  a  may  have  a 
range  or  “window  size”  w ,  such  that  1  <  w  <  r.  With  a  window  size  w,  we  may  have  w 
implicants  that  covers  a.  That  is,  the  minterm  a  can  be  at  any  position  of  the  window. 
For  a  given  axis,  counting  all  window  sizes,  we  have  1+2  +  3  +  ..  .  +  r  =  possible 
implicants  that  cover  a.  Over  the  entire  n-dimensional  space,  we  have  (r^1l)n  =  0(r2n) 
[Ref.  12]. 

B.  THE  PROPERTIES  OF  TRUNCATED  SUM 

There  are  two  important  properties  of  the  truncated  sum  which  are  useful  later  in 
developing  the  ND  algorithm. 

1.  Saturated  minterm  can  be  generated  by  TSUM  operation. 

The  truncated  sum  of  two  or  more  minterms  may  produce  a  saturated 
minterm.  By  definition  4,  the  truncated  sum  of  any  minterm  and  a  sat¬ 
urated  minterm  remains  a  saturated  minterm.  In  other  words,  given  two 
minterms  a,  0  and  minterm  7  =  TSUM(a,0).  If  value  of  7  is  r  -  1,  i.e., 

7  is  a  saturated  minterm  then  for  any  other  minterm  6,  7  +  6  =  7. 
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As  an  example,  in  a  2-variable  4- valued  function,  three  minterms  add 
in  one  position. 


2  xx{  2Xj  +  2  xx}  2x2  4-  1  xxj  2x\  =  3  xxJ  2x\  +  1  1x{  2x2  =  3  xxj  2x 


The  first  two  terms  form  a  saturated  minterm,  and  this  saturated  minterm 
absorbs  the  third  minterm. 

2.  Don't  care  minterm  can  be  produced  by  saturated  minterm. 

In  the  minimization  procedure,  we  may  update  a  minterm  a  to  a' 
by  subtracting  minterm  7  (a'  =  a  —  7),  where  7  is  the  value  of  selected 
implicant.  If  a  G  SAT,  in  succession  of  updates,  the  value  of  a'  may  reach 
the  value  0.  In  that  case,  the  algorithm  will  reset  that  minterm  coordinate 
to  don’t-care ,  i.e.,  value  r. 


C.  DEFINITIONS  USED  IN  ND  ALGORITHM 
Definition  9:  Direct  Neighbors 

Let  a.  and  /?  be  minterms  with  coordinates  <  ai,  a2, . . .  an  >  and  <  bx,  b2, . . .  bn  > 
respectively.  If  for  all  i  we  have  a ,  =  b,  except  one  position  j  such  that  \a3  —  bj\  =  1  we 
say  that  a  and  ft  are  direct  neighbors.  Given  a  minterm  a,  we  use  N(a)  to  denote  the 
set  of  its  direct  neighbors. 

Observation  1:  The  maximum,  number  of  direct  neighbors  of  a  given  minterm  is  2 n. 

Definiton  10:  Directional  Neighbors 

Two  minterms  a  and  axe  directional  neighbors  in  the  direction  xj,  if  a,  =  b,  for 
all  *  €  [l,n]  such  that  i  ^  j.  When  bj  >  a j  we  say  that  (3  is  in  the  positive  direction  of 
a,  while  when  bj  <  aj  we  say  that  0  is  in  the  negative  direction  of  a. 
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Figure  2.2:  Example  for  Connected  Minterms 

Observation  2:  If  0  is  a  direct  neighbor  of  a  then  0  is  a  directional  neighbor  of  a  in 
the  direction  of  Xi  for  some  i  €  [l,nj. 

Definition  11:  Connected  Minterms 

This  is  a  recursive  definition.  Given  a  minterm  a  and  a  minterm  0,  then  we  say  (3 
is  a  connected  minterm  of  a ,  if 

1.  (3  is  a  direct  neighbor  of  a  and  either  g(0)  <  g(a)  or  a  E  SAT. 

2.  0  is  a  directional  neighbor  of  a  in  direction  x,  and  0's  direct  neighbor  is 
connected  to  a  and  either  g(0)  <  g(a)  or  a  €  SAT. 

For  example,  in  Figure  2.2  minterm  2  2x^°x°,  1  3x31X2,  1  2x22X2  and  2  2x23X2 
(pointed  by  arrow)  are  connected  minterms  of  2  2x2  xxJ  (the  minterm  with  @  sign). 

Definition  12:  Connected  Minterm  Count 

CMC  a  is  the  connected  minterm  count  of  minterm  a.  It  is  the  number  of  minterms 
that  are  connected  to  minterm  a. 

Definition  13:  Expandable  Directional  Count 

EDCa  is  the  expandable  directional  count  of  minterm  a.  It  th*»  number  of 
directions  (both  positive  and  negative  for  each  x,)  in  which  a  has  one  or  more  connected 
minterms. 
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Figure  2.3:  Map  for  Example  3,  4;  Step  2  of  Table  3.2 
Observation  3:  0  <  EDCa  <  2n. 

Definition  14:  Clustering  Factor 

The  clustering  factor  relative  to  a  minterm  a  is  defined  as 

CFa  =  EDCa{r-l)  +  CMCa.  (2.4) 

This  is  a  measure  of  the  weight  of  all  connected  minterms  relative  to  a.  The  (r  —  1) 
factor  is  the  range,  or  maximum  possible  number  of  minterms,  in  a  direction  x;. 

Example  3: 

In  Figure  2.1  the  minterm  1  °x°2X2  (the  minterm  with  *  sign)  has  no  connected 
minterms  nor  expandable  directional  neighbors,  i.e.,  its  CMC  and  EDC  values  are  0. 
Figure  2.3  shows  that  a  circled  implicant  1  °xj  1xj  was  subtracted  from  Figure  2.1.  We 
mark  a  minterm  with  a  dot  in  figure  because  it  was  a  saturated  minterm  in  original 
function  map  (see  Definition  8  and  Figure  2.1).  The  minterm  a  =  2lx\  °x°  (the  minterm 
with  @  sign)  is  one  of  the  ten  minterms  and  CMCa  —  3,  EDCa  =  2.  The  clustering 
factors  of  all  minterms  in  Figure  2.3  are  listed  in  Table  2.1. 
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TABLE  2.1:  CFs  for  all  minterms  in  Figure  2.3 


Minterm 

wmaam 

Kimirai 

CF 

9 

4 

6 

18 

18 

Minterm 

gimgga 

CF 

6  1 

13 

4 

10 

10 
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III.  NEIGHBORHOOD  DECOUPLING 

ALGORITHM 


The  neighborhood  decoupling  (ND)  algorithm  [see  Appendix  A  for  the  C  program 
listing]  is  an  improvement  to  Dueck  and  Miller’s  method  [Ref.  8].  Like  the  algorithms 
mentioned  in  Chapter  I  Section  B  and  throughout,  the  ND  algorithm  has  two  compu¬ 
tational  phases:  select  a  minterm  and  select  an  implicant.  Firstly,  the  most  isolated 
minterm  is  chosen  by  using  the  algorithm  M  to  be  described  in  Section  A.  The  most  iso¬ 
lated  minterms  in  general  are  different  from  Dueck  and  Miller’s  method  due  to  different 
decision  rules.  Secondly,  from  all  implicants  which  cover  the  most  isolated  minterm  ND 
algorithm  chooses  the  one  that  is  not  strongly  “coupled”  with  its  neighbors.  In  other 
words,  it  is  the  most  isolated  implicant.  This  decoupling  process  is  based  on  the  ob¬ 
servation  that  if  we  choose  that  specific  implicant  then  we  may  minimize  the  negative 
impact  for  future  minterm  selections  as  well  as  implicant  selections.  The  explanation  of 
this  idea  is  described  in  Section  B.  In  the  algorithm  below,  /  denotes  the  function  to  be 
minimized. 


{ 

SS  <—  <j>;  /*  SS  =  Solution  Set  */ 

WS  —  MSj  =  {a|a  is  generated  by  the  function  /;  if  a  €  SAT  then  mark  its 
coordinate  }. 

While  WS  ±  <j>  do  { 

1.  Use  algorithm  M  (see  Section  A)  to  select  a  minterm  a  from  the  WS. 

2.  Use  algorithm  N  (see  Section  B)  to  select  an  implicant  Ia  that  covers  a. 

3.  SS  «—  SS  U  Ia. 
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4.  V/5  £  Ia  do  { 

compute  g(P)  <-  5(/?)  -  g(a). 
subtracted  Ia  from  W S. 

if  /3  is  originally  marked  and  g{) 3)  =  0  then  g{l 3)  <—  r. 
/*  don’t-care  terms  */ 


} 


} 


} 

A.  ALGORITHM  M:  MINTERM  SELECTION 

Based  on  Definition  15,  all  minterm’s  clustering  factor  can  be  computed.  The  ND 
algorithm  computes  the  clustering  factor  in  the  order  of  coordinates  (x,),  i.e.,  row-column 
order.  For  example,  the  minterm  2  2x\  °x°  is  evaluated  earlier  than  the  minterm  2  2x\  2x\. 
The  algorithm  M  is  described  as  follows: 

1.  Compute  the  corresponding  CFak  for  all  a*  £  WS. 

2.  Select  the  minterm  a  that  has  the  smallest  clustering  factor  (excluding  the  minterm 
a  £  SAT  as  well  as  greater  than  2  x  (r  —  1)).  If  there  is  a  tie,  the  first  one  that 
gets  evaluated  is  chosen. 


B.  ALGORITHM  N:  NEIGHBORHOOD  RELATIVE  COUNT 

The  purpose  of  algorithm  N  is  to  choose  the  most  “isolated”  implicant  ( Ia )  and 
update  the  working  set  IF S'.  It  computes  the  neighborhood  relative  count  (NRC)  for  all 
implicants  that  cover  the  minterm  a.  The  implicant  with  the  smallest  N RC  is  chosen.  In 
other  words,  NRC  is  a  measure  of  the  coupling  strength  of  an  implicant  with  its  neigh¬ 
bors.  To  select  an  implicant  which  is  equivalent  to  breaking  the  coupling  between  that 
implicant  with  its  neighbors,  the  candidate  implicant  should  have  the  smallest  coupling 
strength  with  its  neighbors.  Therefore,  the  ND  algorithm  tends  to  choose  the  most  “iso¬ 
lated”  implicant.  If  there  is  a  tie  in  selecting  the  7a,  ND  algorithm  chooses  the  one  which 
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covers  the  largest  area.  The  computation  of  NRC  for  a  given  implicant  is  described  as 
follows: 

1.  Initialize  the  N RC  to  zero. 

2.  Check  all  neighboring  minterms  of  the  implicant  and  increment  or  decrement  its 
NRC  according  to  the  following  (intuitively  stated)  rule,  which  is,  if  the  coupling 
strength  between  covered  and  uncovered  area  is  weak  (good  for  further  decoupling), 
algorithm  N  decreases  NRC,  otherwise  increases  NRC. 

**********************************************************  ************ 

a:  the  chosen  minterm  from  algorithm  M. 

M:  the  set  of  minterms  which  was  covered  (generated)  by  the  chosen  impli¬ 
cant  ( I  a ). 

N(/3):  the  set  of  direct  neighbors  of  minterm  (3. 

***********************************************************************  * i 

{ 

NRC  -  0; 

V,3  €  M  and  (3  ^  a  do  { 

if (g(P)  -  g(a)  <  0)  then  NRC  *-  NRC  -  2; 

} 

V/?  €  M  and  V7  €  N(/3)  do  { 

if  (7  &  M  and  7  ^  0  and  (7  £  SAT  or  (3  £  SAT))  then  { 
if  W)  ~  ff(o)  >  9(7)  )  then  { 

if  (7  €  SAT  )  then  NRC  «-  NRC  -  1; 

else  NRC  «-  NRC  +  2; 

} 

if  W)  ~  9(<x)  <  9(lf)  )  then  { 

if  {g(P)  -  9{ 7)  )  then  NRC  <-  NRC  +  2; 
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if  (7  €  SAT  and  5(7)  -  g(3 )  <  0  )  then 
NRC  *-  NRC  +  2; 

else  { 

if  (g(/3)  >  g(a)  and  g(/3)  ^  g{ 7)  )  then  { 
if(  3  €  SAT )  then  NRC  <-  NRC  -  1; 

else  NRC  *-  NRC  +  2; 

}  /*  end  if  */ 

}  /*  end  else  */ 

}  /*  end  if  7 


if  ( g(/3 )  -  s(a)  =  5(7)  )  then  { 

if  (<7(7)  >  0  or  f3  €  SAT)  then 

JVflC  NRC  -  1; 
else  NRC  iV/?C  -  2; 


} 

}  /*end  if  7 

} 

if  (A/  =  {o})  then  { 

if  (a  €  SAT)  then  NRC  <—  2; 

else  if  ( NRC  <  0)  then  NRC  *-  1; 

} 

else  NRC  NRC  +  2; 
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Figure  3.1:  Map  for  Example  4;  Step  3  of  Table  3.2 

Example  4: 


It  is  instructive  to  examine  an  example  of  the  application  of  the  ND  algorithm. 
We  use  a  two-variable  four-valued  function  as  an  example.  The  input  function  to  be 
minimized  is  expressed  as: 


f  _  q  0  3  11  ,  o  1_2  0  0  I  q  11  2  3  ,  o  2„,2  2  3  ,  i  2^2  3  3  ,  .  0  0  2  2 
J  —  o  X |  X2  ~T  £  Xj  X2  1  <3  Xj  X2  t  ^  Xj  X2  I  X^  X2  “I"  1  Xj  X2» 


The  working  set,  VT5,  is  initialized  to  MSf  and  is  represented  in  Figure  2.1.  The  clus¬ 
tering  factors  of  all  minterms  in  VF5  are  calculated  (see  example  3  for  computation). 
The  smallest  CF  that  first  comes  from  minterm  1  °Xi  2x\  ;  therefore,  algorithm  M  will 
select  a  =  1  °Xj2X2.  The  ND  algorithm  computes  the  NRCi  for  each  implicant  /  which 
covers  a  using  algorithm  N.  Since  implicant  1  °xj  'x]  has  the  smallest  NRC(— 2),  the 
ND  algorithm  selects  it  as  the  first  implicant  in  the  solution  set  ( SS ).  Table  3.2,  together 
with  Figures  2.3,  3.1,  and  3.2  shows  the  steps  of  choosing  successive  implicants.  The 
*  sign  in  each  figure  indicates  the  most  isolated  minterm  while  a  circled  implicant  is 
the  most  isolated  implicant.  Suppose  we  have  chosen  two  implicants  from  the  function 
shown  in  Figure  2.1  the  resulting  function  is  shown  in  Figure  3.1.  The  minterm  1  2x\ 3x\ 
is  selected,  since  it  has  the  smallest  CF{ 4). 
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Figure  3.2:  Map  for  Example  4;  Step  4  of  Table  3.2 

TABLE  3.1:  NRCs  for  Minterm  20x°xx\  in  Figure  3.2 


Implicant 

2  V  W 

NRC 

2 

0 

-2 

-4 

There  are  two  implicants  that  cover  the  minterm  1  2x\  3x3,  and  their  NRC  values 
are  2  and  —2.  The  implicant  1  1ij3x3  is  chosen  since  it  has  the  smallest  NRC(— 2). 
Having  updated  the  working  set  and  added  1  xx\*x\  to  the  solution  set,  we  have  the 
new  map  in  Figure  3.2.  From  Figure  3.2,  the  NRC  values  for  the  minterm  2  °x°  'xj  are 
available  in  Table  3.1.  Finally,  the  working  set  should  contain  value  0  (empty  square)  or 
4  ( don’t-care )  as  shown  in  Figure  3.3. 


Figure  3.3:  Final  Working  Set 
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TABLE  3.2:  Steps  of  ND  algorithm 


Step 

Minimum 

CFa 

Minterm  a 

Minimum 

NRC 

Candidate  Implicant 

Reference 

1 

0 

-2 

2  °x11 1x22 

Figure  2.1 

2 

4 

reaami 

-8 

Figure  2.3 

3 

4 

mvmmm 

-2 

1  1x153x23 

Figure  3.1 

4 

6 

-4 

3VV 

Figure  3.2 

The  final  minimized  result,  g,  is  expressed  as: 


9 


=  2  °xj  xx\  +  1 


Jx2  +  1  xx\zx\  +  3 


0 


X 


3  1 
1 


X 


1 

2' 
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IV.  COMPARISON  RESULTS 

In  this  thesis  all  testing  results  were  obtained  by  running  sample  functions  on  the 
VAX  11/785  and  ISI  workstations.  A  large  number  of  sample  functions  were  randomly 
generated  (63,500).  We  applied  each  algorithm  to  minimize  these  sample  functions  in 
a  way  similar  to  the  method  used  by  Tirumalai  and  Butler  [Ref.  9],  and  Yurchak  and 
Butler  [Ref.  10].  However,  Tirumalai  and  Butler  generated  the  input  functions  with 
a  fixed  number  of  minterms.  We  generated  the  input  function  with  a  fixed  number  of 
product  terms.  This  thesis  investigated  three  algorithms:  (1)  Pomper  and  Armstrong 
[Ref.  5],  (2)  Dueck  and  Miller  [Ref.  7,  8,  10],  and  (3)  Neighborhood  Decoupling  under 
various  settings: 

1.  Two- variable  four- valued  with  3  to  16  input  product  terms. 

2.  Two-variable  five-valued  with  3  to  25  input  product  terms. 

3.  Three-variable  four-valued  with  3  to  45  input  product  terms. 

4.  Four-variable  four-valued  with  3  to  35  input  product  terms. 

For  two-variable  four-valued  functions,  we  generated  1000  sample  functions  for  a 
given  number  of  input  product  terms  from  3  to  16,  i.e.,  we  generated  14,000  functions 
to  test  the  different  algorithms.  For  the  other  three,  we  generated  500  functions  for  each 
given  number  of  input  product  terms,  i.e.,  for  each  set  we  generated  11,500,  21,500,  and 
16,500  functions  respectively. 

In  Section  B,  functions  with  a  different  number  of  product  terms  are  tested  and  the 
computation  times  are  recorded.  For  a  given  number  of  product  terms,  the  computation 
time  from  two  to  four  variable  functions  grows  exponentally.  This  explains  why  we  did 
not  produce  more  than  1000  functions  for  each  case. 
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A.  PERFORMANCE  COMPARISON 


There  are  two  ways  to  measure  the  performance.  First,  the  average  number  of 
output  product  terms  will  show  us  the  advantage  of  each  algorithm.  Second,  the  ratio 
of  the  performance  results  will  show  us  the  performance  increases  or  decreases  for  each 
algorithm  under  different  conditions. 

1.  AVERAGE  NUMBER  OF  OUTPUT  PRODUCT  TERMS 
For  each  set,  we  compute  the  average  number  of  output  product  terms  (see  Ap¬ 
pendix  B).  From  these  data  points,  we  can  plot  a  curve,  indicating  the  average  number 
of  output  product  terms  as  a  function  of  input  product  terms.  The  plots  for  the  various 
sets  are  shown  in  Figures  4.1,  4.2,  4.3,  and  4.4.  In  these  bell-shaped  figures  we  observed 
three  important  features: 

1.  The  differences  in  the  number  of  average  output  product  terms  among  those  algo¬ 
rithms  axe  insignificant,  when  the  input  number  of  product  terms  are  small.  When 
the  number  of  input  product  terms  is  small,  the  function  is  simple  and  easy  to 
minimized;  therefore,  the  number  of  output  product  terms  for  each  algorithm  are 
almost  same. 

2.  The  highest  point  of  each  bell-shaped  curve  is  the  highest  average  number  of  output 
product  terms  for  each  algorithm.  The  number  of  input  product  terms  that  has  the 
highest  average  number  of  output  product  terms  constitutes  the  most  complicated 
functions. 

3.  When  the  number  of  input  product  terms  get  larger,  the  curve  goes  down  gradually. 
This  is  chiefly  due  to  the  fact  that  the  more  input  product  terms  in  a  function  the 
higher  tendency  of  generating  saturated  minterms.  In  most  cases,  a  single  implicant 
can  cover  a  cluster  of  saturated  minterms. 


From  Figures  4.1,  4.2,  4.3,  and  4.4,  the  neighborhood  decoupling  algorithm  obvi¬ 
ously  outperforms  the  other  two  algorithms. 
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Number  of  Input  Product  Terms 

Figure  4.1:  Two- Variable  Four- Valued  Average  Product  Term 
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Figure  4.3:  Three- Variable  Four- Valued  Average  Product  Term 


5  10  15  20  25  30 

Number  of  Input  Product  Terms 

Figure  4.4:  Four- Variable  Four- Valued  Average  Product  Term 


TABLE  4.1:  Test  Results  of  14000 

2- Variable  4-Val 

!ued  Sample  Functions 

Performance 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood  Decoupling 

54 

212 

336 

479 

2005 

2308 

equal 

11002 

11002 

11002 

total 

11535 

13219 

13646 

ratio 

0.8239 

0.9442 

0.9747 

2.  PERFORMANCE  RATIOS 

We  consider  another  measure,  the  performance  ratio,  to  demonstrate  the  perfor¬ 
mance  of  each  algorithm.  First,  for  each  algorithm,  we  count  the  number  of  minimized 
functions  that  specify  which  algorithm  is  the  “best”  of  the  three,  i.e.,  the  number  of  in¬ 
stances  that  a  specific  algorithm  uses  the  minimum  number  of  implicants  (output  product 
terms).  If  two  algorithms  use  an  equal  number  of  implicants  and  less  than  the  other  one, 
we  determine  that  they  are  “better”  than  the  third  one.  When  all  three  algorithms  use  an 
equal  number  of  implicants  to  minimize  a  function,  we  determine  that  they  are  “equal”. 
The  performance  ratio  is  defined  as 


T]  —  (Nbcst  +  Nbette r  +  ^equal)  4-  Ntotal,  (4-1) 

where  Nb„t ,  N better i  and  Nequai  are  the  number  of  instances  that  specific  algorithm  per¬ 
forms  “best”,  “better”,  and  “equal”  respectively. 

The  total  number  of  functions  tested,  Ntotai  in  our  case  is  14,000,  11,500,  21,500, 
or  16,500.  Tables  4.1,  4.2,  4.3,  and  4.4  show  the  performance  ratios  for  each  setting. 
For  example,  in  Table  4.1,  with  14,000  functions  tested,  we  counted  the  cases  which 
Neighborhood  Decoupling  algorithm  performs  no  worse  than  the  others  as  13,646.  That 
is,  Vn:  —  13646/14000  =  0.9747.  Tables  4.1,  4.2,  4.3,  and  4.4  show  that  the  performance 
is  degraded  when  n  or  r  is  increased.  However,  the  neighborhood  decoupling  algorithm 
obviously  outperforms  the  others. 
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TABLE  4.2:  Test  Results  of  1150C 

2- Variable  5-Va 

ued  Sample  Functions 

Performance 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood  Decoupling 

best 

114 

356 

1248 

better 

762 

1993 

2509 

equal 

7150 

7150 

7150 

total 

8026 

9499 

10907 

ratio 

0.6979 

0.8260 

0.9484 

TABLE  4.3:  Test  Results  of  2150C 

3- Variable  4-Va! 

ued  Sample  Functions 

Performance 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood  Decoupling 

948 

1595 

4578 

2791 

4307 

5918 

equal 

7871 

7871 

7871 

total 

11610 

13773 

18367 

ratio 

0.5400 

0.6406 

0.8543 

TABLE  4.4:  Test  Results  of  1650C 

i  4- Variable  4-Val 

ued  Sample  Functions 

Performance 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood  Decoupling 

1604 

3077 

4952 

1420 

3106 

3440 

equal 

2884 

2884 

2884 

total 

5908 

9067 

11276 

ratio 

0.3581 

0.5495 

0.6834 
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B.  TIMING  COMPARISON 


HAMLET  CAD  (Computer  Aided  Design)  tool  has  an  option  that  can  measure 
program’s  computation  time.  In  this  section,  the  algorithm  timing  comparison  counts  the 
average  computation  time  for  each  input  product  term  minimization.  These  measures 
are  collected  on  the  ISI  workstations  at  NPS.  The  theoretical  worst  case  computional 
complexity  of  each  algorithm  is  0(r2n)  (see  Chapter  II).  The  decision  rules  for  selecting 
minterm  a  and  implicant  Ia  in  each  algorithm  are  different.  In  general,  an  algorithm 
which  needs  more  or  complex  decision  rules  in  selecting  a  and  Ia  should  take  longer 
computation  time.  There  may  be  a  trade-off  between  performance  and  computing  time. 
In  each  computation  cycle  (select  a  minterm  and  select  an  implicant),  each  algorithm 
generates  an  output  product  term.  The  ND  algorithm  uses  more  complex  decision  rules 
than  Dueck  and  Miller’s.  It  may  appear  that  ND  algorithm  should  use  more  computation 
times.  However,  this  is  not  true,  since  ND  algorithm  in  general  has  less  computation  cycle 
than  other  algorithms  so  that  it  stop  computation  earlier  than  other  algorithms.  From 
the  analysis  of  testing  case  results  (see  Appendix  C),  the  average  computation  time  of 
the  ND  algorithm  is  no  worse  than  Dueck  and  Miller’s  for  more  than  4  second,  sometimes 
even  better. 
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Average  Computation  Time  (sec)  Average  Computation  Time  (sec) 


Neighborhood  Decoupling  x; 

Dueck  &  Miller  o 
Pomper  &  Armstrong  * 


Number  of  Input  Product  Terms 

Figure  4.5:  Two- Variable  Four- Valued  Average  Computation  Time 


Number  of  Input  Product  Terms 

Figure  4.8:  Two- Variable  Five- Valued  Average  Computation  Time 
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V.  DISCUSSIONS  AND  CONCLUSIONS 


A.  DISCUSSIONS 

Notice  that  in  the  comparisons  we  did  not  include  the  Bessilich’s  algorithm  and  the 
absolute  minimum  solutions.  We  now  justify  the  reasons.  Besslich’s  algorithm  was  not 
included  in  HAMLET  because  it  requires  a  truth  table  to  operate  efficiently.  We  did  not 
have  room  to  store  large  truth  tables.  Also  Besslich’s  algorithm  is  likely  to  require  much 
time  [Ref.  10].  In  addition,  we  know  that  the  Dueck  and  Miller  is  a  satisfactory  [Ref.  9, 
10]  heuristic  which  is  reasonable  approximation  to  the  absolute  minimum  solutions  for 
two-variable  four-valued  functions. 

The  performance  results  show  that  the  ND  algorithm  outperforms  the  other  two 
algorithms.  Although  Pomper  and  Armstrong’s  algorithm  is  faster  than  the  other  two 
algorithms,  about  one  third  of  the  computation  time,  the  number  of  average  output 
product  terms  and  performance  ratios  are  worse  than  the  other  two  algorithm’s.  From 
the  timing  comparison  results,  the  ND  algorithm  runs  faster  than  Dueck  and  Miller’s  in 
larger  variable  setting  conditions  (more  than  two  variable).  This  time  efficiency  is  due 
to  the  decision  rules  employed  in  the  ND  algorithm  which  take  advantage  of  the  special 
property  of  truncated  sum  operations.  In  other  words,  the  input  product  terms  have  a 
tendency  to  produce  saturated  minterms.  In  the  decoupling  process  (algorithm  N),  a 
minterm  in  SAT  will  always  qualify  to  combine  with  its  neighbors  to  form  an  implicant. 
Although  the  ND  algorithm  conceptually  is  similar  to  Dueck  and  Miller’s  algorithm,  the 
ND  algorithm  uses  saturated  minterms  in  an  more  effective  way.  For  example,  when 
we  update  (deduct)  saturated  minterm  from  expression,  the  minterm  will  be  updated  to 
a  “ don’t-care  minterm”  (see  Chapter  III).  Like  binary  logic  minimization,  a  “ don’t-care 
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minterm”  [Ref.  13]  can  simplify  the  minimization  process. 

Recall  that  in  algorithm  N,  we  compute  the  NRC  values  for  a  given  minterm  a  by 
examining  the  relationships  of  Ia  and  its  immediate  neighbors,  i.e.,  one  step  look-ahead. 
It  is  natural  to  believe  that  with  more  steps  look-ahead  we  might  make  a  better  choice 
of  the  implicant  and  therefore  provide  a  better  solution.  The  exponential  growth  of  the 
number  of  possible  implicants  restricts  the  practical  use  of  k-lookahead  for  k  1. 

B.  CONCLUSION 

The  truncated  sum  MVL  minimization  can  be  done  by  the  neighborhood  decou¬ 
pling  algorithm  which  selects  the  most  isolated  minterms  as  well  as  implicants.  In  the 
development  of  ND  algorithm,  we  have  the  following  observations:  (1)  truncated  sum 
operations  may  produce  saturated  minterms,  and  (2)  a  saturated  minterm  will  reduce  to 
a  don’t-care  minterm  in  the  minimization  process.  The  decision  rules  of  the  ND  algorithm 
take  full  advantage  of  these  observations.  The  ND  algorithm  outperforms  most  heuristic 
methods  and  does  not  lose  its  run  time  efficiency  because  the  algorithm  finds  the  solution 
and  stops  earlier  than  others.  As  the  number  of  variable  of  the  input  function  become 
larger,  i.e.,  greater  than  three,  the  ability  to  minimize  the  function  in  reasonable  time  is 
an  important  and  challenging  research  area. 
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APPENDIX  A:  PROGRAM  LISTINGS 


/*  nd.c 

-  This  module  implements  the  Neighborhood  Decoupling  heuristic  for 
MVL  minimization. 

-  This  program  can  be  added  to  HAMLET  [Ref .  10] . 

-  HAMLET  is  written  by  the  instructor  J.  Yurchak  in  the  Department 
of  Computer  Scinece  at  Naval  Postgratuate  School. 

*/ 

•include  "defs.h" 
static  int 

better.found; 


/*  Neighborhood  Decoupling  functions  -  */ 

N_D() 

/* . - . - . 

: function: 

-  Perform  the  Neighborhood  Decoupling  heuristic  on  the  input  expression 
: algorithm: 

Start  with  a  working  copy  E.work  of  the  original  function  E_orig; 
Initialize  a  final  function  E.final; 

While  (there  are  still  minterms  to  pick)  { 

Pick  a  minterm  X  from  E.work; 

Pick  the  best  implicant  I  for  X; 

Add  I  to  E.final; 

Subtract  I  from  E.work; 

> 

:globals: 

E.orig 

e-flag 
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m_flag 
q-flag 
G_f lag 
FO.ratio 
:side_ef facts: 

STAT 
HEUR 
E_work 
E_f  inalP 
:called_by : 

mainO 
: calls: 

daalloc.exprO 

dup.exprO 

print_terms() 

print_map() 

mim() 

pick_implicant() 

subtract_implicant() 

print_source() 

. - . - . */ 

int  num_ impl  =  0; 

int  *X ; 

Imp li cant  *1; 

float  ratio; 

if  (E_final[ND].I  !=  NULL) 

dealloc_expr(4E_f inal [ND] ) ; 

fifdef  ANALYZER 

STAT  *  ftWY.stat; 

fendif 

HEUR  *  ND; 

dup.expr (AE.work , AE.orig) ; 

E_f inalCHEUR] .nterm  ■  0; 

E_f inalCHEUR] .radix  ■  E.orig. radix; 

E_f inal[HEUR] .nvar  *  E.orig.nvar; 

E_f inalCHEUR] .1  «  NULL; 

fifdef  ANALYZER 
if  (e„flag) 

print_terms(tE_orig) ; 
if  (m.flag)  { 
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printf ("  Orig  map  (nda):\n"); 
print _map() ; 

> 

tendif 


better.found  =  Sc.flag; 
resource.used(START) ; 
for  (; ;)  { 

if  ( (X  =  mim(tE.work))  ==  NULL)  { 
if  (num_impl  <  E_orig.nterm) 
better.found  =  1; 
break; 

> 

I  =  pick_implicant(X) ; 
num_impl++; 

subtract. implicant (I) ; 

#ifdef  ANALYZER 

if  (i.flag) 
print_implicant(X,I) ; 
if  (m.flag) 
print.mapO ; 

tendif 

if  (Sm.flag)  { 

if  (num.impl  >=  E.orig .nterm) 
break; 

> 

> 

resource_used(STOP) ; 

if  ( !better_found)  { 

num.impl  *  E.orig.nterm; 
dup_expr(A(E_f inal[ND] ) ,*E_orig) ; 


iifdef  ANALYZER 

if  ( Iq.flag  kt  'G.flag)  { 
if  (! better.found) 

printf ("  X-4d  nda:  X4d  /X-4d  X61d.7.2.21d\n\n", 

expr.seq, num.impl, num.impl ,secs_used() ,tsecs_used() ) ; 

else 


printf ("  X-4d  nda:  X4d  /X-4d  X6d.X2.21d\n\n" , 


expr_seq,num_impl,E_orig.nterm, 
secs.usedO ,tsecs_used()) ; 

> 

tendif 


dealloc_expr(4E_work) ; 

> 


static  in  *mim(E) 

Expression  *E; 

/* . - . - — . 

: function: 

-  Compute  the  clustering  factor  and  find  the  Most  Isolated  Minterm 
in  the  expression  pointed  to  by  E,  and  return  its  coordinates  as 
a  vector. 

-  Local  to  nd.c 
••globals: 

radix 

nvar 

:side_effects: 

STAT 

:called_by: 

N.DO 

: calls: 

next_coord() 
evalO 
vcopyO 
: returns : 

-  A  vector  of  integers  representing  the  coordinate  of  the  most 
isolated  minterm,  or  NULL  if  no  more  minterms. 

-  The  value  at  that  location  is  also  returned  as  the  last  integer 
in  the  vector. 

. - . . . - . */ 

{ 

register  i,j,k; 

int  cur.val  ■  E->radix, 

cur.CF  ■  MAX. INT, 

X_orig[MAX_VAR+2] , 

R_1  *  radix  -  1, 

Not.all  ■  0, 

All _t run  ■  0, 

TRUN  »  2*R.l, 
last  ■  0, 
expanded, 
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value [2] , 

vail [2] , 

val2 [2] , 

cf, 

ea, 

dea, 

term; 

int  *X,*next_coord() ; 

static  int 

coord[MAX_VAR+2]  , 
save_coord[MAX_VAR+2] ; 

#if def  ANALYZER 

STAT->calls_mim++ ; 

#endif 

for  (term=0;  term  <  E_orig.nterm;  term++)  { 
k  =  1; 

while  ((X=next_coord(coord,ft(E->I[term] ) ,k))  !=  NULL)  { 
vcopy(value,eval(E,X)) ; 

if  (value [EVAL]  kk  value [EVAL]  <  radix)  { 
if  (! value [HLV]) 

Not.all  *  1; 
if  (All.trun)  { 
cf  «  0; 
dea  *  0; 
ea  =  0; 

for  (j=0 ;  j  <  nvar;  j++)  X_orig[j]  =  X[j]  ; 
/*  for  each  variable  (direction) . . .  */ 
for  (j*0;  j  <  nvar;  j++  )  { 
expanded  »  0; 

/*  If  not  on  a  left  hand  edge,  move  left  */ 
while  (X [j]  >  0)  { 

xCj]~; 

vcopy(vall,eval(E,X)) ; 

if  (vail [EVAL])  { 
expanded  m  1; 
ea++; 

> 

else 

break; 

> 

X[j]  ■  X.origCj]; 

if  (expanded)  { 
expanded  *  0; 
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dea++; 


> 

/*  if  we  didn’t  start  on  a  right  hand  edge,  move  right  */ 
while  (X[j]  <  R.l)  { 

X[j]++; 

vcopy(val2,eval(E,X)) ; 
if  Cval2[EVAL] )  { 
expanded  =  1; 
ea++; 

> 

else 

break ; 

> 

X[j]  =  X.origCj]; 
if  (expanded) 
dea++ ; 

> 

/*  compute  the  clustering  factor  */ 
cf  =  (dea  *  R.l)  +  ea; 
if  (cf  <  cur.CF)  { 

cur.val  *  value [EVAL] ; 
cur.CF  *  cf; 

for  (i*0;  i  <  nvar;  i++)  save_coord[i]  •  X[i] ; 

> 

> 

else  { 

cf  *  0; 
dea  *  0; 
ea  *  0; 

for  (j=0;  j  <  nvar;  j++)  X.origCj]  =  X[j] ; 

/*  for  each  variable  (direction) . . .  */ 
for  (j=0;  j  <  nvar;  j++  )  •( 
expanded  ■  0; 

/*  If  not  on  a  left  hand  edge,  move  left  */ 
while  (X[j]  >  0)  { 

X[j]~; 

▼copy(vali,eval(E,X)) ; 

if  (vail [EVAL]  At  (vail [EVAL]  <=  value [EVAL] 
II  value [HLV]))  { 
expanded  •  1; 
ea++; 

> 

else 
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break; 

> 

XCj]  =  X.origCj]  ; 
if  (expanded)  { 
expanded  =  0; 
dea++; 

> 

/*  if  we  didn’t  start  on  a  right  hand  edge,  move  right  */ 
while  (X[j]  <  R_l)  { 

X[j]++; 

vcopy(val2,eval(E,X)) ; 

if  Cval2[EVAL]  44  (val2[EVAL]  <=  value [EVAL] 
I  I  value [HLV] ) )  { 
expanded  =  1; 
ea++; 

> 

else 

break; 

> 

XCj]  =  X.origCj]; 
if  (expanded) 
dea++; 

> 

/*  compute  the  clustering  factor  */ 
cf  *  (dea  *  R.l)  +  ea; 

if  ( ! (value [HLV]  44  cf  >  TRUN))  { 
if  (cf  <  cur.CF)  { 

cur.val  =  value [EVAL] ; 
cur.CF  »  cf; 

for  (i*0;  i  <  nvar;  i++)  save_coord[i]  *  X [i] ; 

> 

> 

> 

>  /*  end  if  */ 

k  *  0; 

>  l*  end  while  */ 

if  (Hast  44  (term  «■  (E.orig.nterm  -  1))  44  INot.all)  { 

All.trun  ■  1; 
cur.CF  ■  MAX.INT; 
term  ■  -1; 
last  *  1; 
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if  (cur.CF  ==  KAX.IHT) 
return(NUIL) ; 

save_coord[nvar+l]  =  cur.CF; 
save_coord[nvar]  =  cur.val; 

return(save_coord) ; 


static  int  valid_implicant(I) 

Implicant  *1; 

/* . - . 

: function: 

-  Decide  upon  the  validity  of  implicant  I 

-  Local  to  nd.c 

:globals : 

E.work 

E_orig 

:side_effects: 

STAT 

: called _by : 

pick_implicant() 

: calls: 

next_coord() 

eval() 

vcopyO 

: returns : 

1  if  a  valid  implicant 
0  if  not 


int 

*X, 

init  *  1, 

R_1  •  radix  -  1, 
value  ■  I->coeff, 

Vo  [2]  ,Vw[2]  ; 
static  int 

coord [MAX. VAR+2] ; 

fifdef  ANALYZER 

STAT->calls_valid_implicant++ ; 

iendif 
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while  ((X  =  next_coord(coord,I,init))  !=  NULL)  { 
init  =0; 

vcopy  (Vw ,  eval  ( StE.work ,  X) ) ; 
vcopy(Vo,eval(4E_orig,X)) ; 

if  (((Vw [EVAL]  <  value)  kk  «Vw[HLV])  kk  (Vo [EVAL]  <  R.l)) 
return(O) ; 

> 

return(l) ; 

> 


static  int  compute_nrc(I) 

Implicant  *1; 

/* . . . — 

.•function: 

-  Compute  the  NRC  for  the  given  implicant 

-  Local  to  nd.c 
:globals : 

radix 

nvar 

: side.eff ects : 

STAT 

:called_by: 

pick. implicant ( ) 

: calls: 

next.coordO 
evalO 
vcopy () 

: returns : 

-  an  integer  NRC 


{ 

int  *X ; 

int  I.value  *  I->coeff; 

register  i; 

int  value [2] , 

R.l  *  radix  -  1, 

neighbor. value [2] , 

good, 

bad, 

diff , 

equal , 

neig.boun, 

first, 

nrc  *  0, 


*/ 
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init  *  1; 
static  int 

coord [MAX.VAR+2] ; 

#ifdef  ANALYZER 

STAT->calls_compute_nrc++ ; 

# end if 


/*  for  each  coordinate  in  the  jmplicant  . . .  */ 
while  ((X  =  next_coord(coord,I,init))  !=  NULL)  { 
init  =  0; 
equal  =  0; 

vcopy(value,eval(4E_work,X)) ; 
if  (value [EVAL]  *=  radix) 
continue; 

diff  =  value [EVAL]  -  I.value; 
first  *  1; 

/*  for  each  direction  ...  */ 

for  (i=0;  i  <  nvar;  i++)  { 
good  ■  0; 
bad  *  0; 

if  ((diff  <=  0)  kk  first)  { 
good  *  2 ; 

first  *  0; 

> 

/*  if  there  is  a  left  neighbor,  examine  it  */ 

if  (X[i]  !=  0  kk  X[i]  **  I->B[i] .lower)  { 

X[i]  — ; 

vcopy (neighbor. value , eval (ftE.work , X) ) ; 
neig.boun  *  neighbor. value [EVAL]  -  value [EVAL]; 
X[i]++; 

if  (neighbor .value [EVAL]  !=  0)  { 

if  (!neighbor_value[HLV]  II  ! value [HLV] )  { 
if  (neighbor_value[EVAL]  <  diff)  { 
if  (neighbor_value[KLV] ) 
good  +*  1; 

else 

bad  +■  2; 

> 

if  (neighbor.value [EVAL]  >  diff)  { 
if  (! neig.boun) 
bad  +»  2; 

if  (neighbor.value [HLV]  kk  neig.boun  <  0) 
bad  +■  2; 

if  (diff  >  0  kk  neig.boun)  { 
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if  (value [HLV]) 
good  +=  1 ; 

else 

bad  +=  2; 

> 

> 

else  -C 

if  (neighbor_value[HLV]  ||  value  [HLV]) 
good  +=  1; 

else 

good  +=  2; 

> 

> 

>  /*  end  if  */ 

> 

/*  if  there  is  a  right  neighbor,  examine  it  */ 

if  (X[i]  !=  R_1  kk  X[i]  --  I->B[i] .upper)  { 

X[i]++; 

vcopy (neighbor. value , aval (ftE.work , X) ) ; 
neig.boun  ■  neighbor.value [EVAL]  -  value [EVAL] ; 

X  [i]  — ; 

if  (neighbor.value [EVAL]  !=  0)  { 

if  (! neighbor. value  [HLV]  ||  ! value [HLV])  { 
if  (neighbor.value [EVAL]  <  diff)  { 
if  (neighbor.value [HLV] ) 
good  +=  1; 

else 

bad  +=  2; 

> 

if  (neighbor.value [EVAL]  >  diff)  { 
if  (Ineig.boun) 
bad  +=  2; 

if  (neighbor.value  [HLV]  kk  neig.boun  <  0) 
bad  +*  2; 

if  (diff  >  0  kk  neig.boun)  { 
if  (value [HLV]) 
good  +■  1; 

else 

bad  +*  2; 

> 

> 

else  { 

if  (neighbor.value [HLV]  ||  value [HLV]) 
good  ♦»  1; 


else 


> 


> 


good  +=  2; 


/*  update  the  nrc  */ 

nxc  *  (nrc  -  good)  +  bad; 
>  /*  end  for  */ 

>  /*  end  while  */ 
return (nrc) ; 


static  Implicant  *pick_implicant(X) 

int  *X ; 

/*  . - . . 

: function: 

-  Pick  the  best  implicant  for  minterm  X 

:globals : 

radix 

:side_effects: 

STAT 

:called_by: 

N_D() 

: calls: 

init_ implicant () 
gen .bounds () 
next _ impl icant ( ) 
evalO 
vcopy() 
compue.nrc ( ) 
copy.implicantO 
valid.implicantO 

: returns : 

-  A  pointer  to  a  term  representing  the  best  implicant. 


{ 

int  cur.nrc  *  MAX.INT, 

nrc  ■  0, 

I .value, 
i. 

init  ■  1, 
first  »  i; 


*/ 
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Implicant 


*i; 

static  int 
coord [MAX. VAR+2]  ; 
static  Bound 

I.bound [MAX. VAR+2]  ; 
static  Implicant 
I.best; 

Bound  *B ; 

int  V  [2] , 

value [2] ; 

iifdef  ANALYZER 

STAT->calls_pick_implicant++ ; 

#endif 

I.best. B  =  I.bound; 
init.implicant(X) ; 

B  *  gen.bounds(X) ; 
vcopy(V,eval(AE_orig,X)) ; 
while  ((I  *  next.implicant(B))  !*  NULL)  { 
if  (V  [HLV]  )  { 

for  (I->coeff*X[nvar] ;  I->coeff  <  radix;  (I->coeff )++)  { 
if  (valid.implicant(I))  { 
nrc  ■  compute.nrc(I) ; 
if  (first) 
nrc  =  2; 

else 

nrc  +=  2; 

if  (nrc  <=  cur.nrc)  { 
cur.nrc  3  nrc; 

I->nrc  =  nrc; 

copy. implicant (tl.best ,1) ; 

> 

> 

> 

first  ■  0; 

> 

else  { 

I->coeff  ■  X[nvar]; 
if  (valid.implicant(I))  { 
nrc  ■  compute.nrc(I) ; 
if  (first)  { 
first  »  0; 
if  (nrc  <  0  ) 
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nrc  =  1; 


nrc  +=  2; 


nrc  +=  2; 

if  (nrc  <=  cur.nrc)  { 
cur.nrc  =  nrc; 

I->nrc  =  nrc; 

copy_implicant(4I_best,I) ; 

> 


return (ftl .best) ; 
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APPENDIX  B:  AVERAGE  NUMBER  OF 
OUTPUT  PRODUCT  TERMS 


TABLE  B.l:  Two- Variable  Four- Valued  Average  Output  Product  Terms 


Number  of 
Input  Terms 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood 

Decoupling 

3 

2.838 

2.716 

2.709 

4 

3.483 

3.290 

3.277 

5 

3.916 

3.710 

3.675 

6 

4.178 

3.959 

3.934 

7 

4.203 

4.023 

3.977 

8 

4.201 

4.009 

3.950 

9 

4.072 

3.920 

3.872 

10 

3.913 

3.763 

3.699 

11 

3.717 

3.583 

3.534 

12 

3.573 

3.438 

3.398 

13 

3.355 

3.233 

3.207 

14 

3.178 

3.079 

3.042 

15 

2.991 

2.899 

2.877 

16 

2.759 

2.696 

2.658 

45 


TABLE  B.2:  Two- Variable  Five- Valued  Average  Output  Product  Terms 


Number  of 
Input  Terms 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood 

Decoupling 

3 

3.122 

2.858 

2.832 

4 

4.044 

3.718 

3.684 

5 

4.660 

4.296 

4.236 

6 

5.222 

4.868 

4.760 

7 

5.540 

5.138 

5.040 

8 

5.740 

5.428 

5.258 

9 

5.832 

5.548 

5.360 

10 

5.948 

5.566 

5.420 

11 

5.814 

5.530 

5.372 

12 

5.762 

5.524 

5.306 

13 

5.486 

5.246 

5.062 

14 

5.364 

5.198 

4.894 

15 

5.014 

4.858 

4.622 

16 

4.780 

4.668 

4.376 

17 

4.484 

4.360 

4.172 

18 

4.250 

4.184 

4.016 

19 

4.032 

3.942 

3.762 

20 

3.822 

3.766 

3.600 

21 

3.730 

3.676 

3.504 

22 

3.564 

3.482 

3.362 

23 

3.224 

3.152 

3.032 

24 

3.044 

2.998 

2.902 

25 

2.982 

2.944 

2.840 
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TABLE  B.3:  Three- Variable  Four- Valued  Average  Output  Product  Terms 


Number  of 
Input  Terms 


3 

4 

5 

6 


3.026 

4.090 

5.162 

6.168 

7.186 

7.916 

8.466 

9.02 

9.646 
10.056 
10.426 
10.588 
10.730 
10.668 
10.676 
10.866 
10.628 
10.586 
10.508 
10.280 
10.164 
9.726 
9.448 
9.142 
8.896 
8.706 
8.820 
8.648 
8.592 
8.104 
8.022 

7.646 
7.160 
7.076 
6.888 
6.632 
6.492 
6.406 
6.128 
5.878 
5.570 
5.270 
5.302 


Dueck  and  Miller 

Neighborhood 

Decoupling 

2.940 

2.928 

3.892 

3.862 

4.834 

4.840 

5.756 

5.718 

6.578 

6.492 

7.372 

7.238 

7.852 

7.762 

8.570 

8.324 

9.068 

8.814 

9.578 

9.266 

9.918 

9.608 

10.126 

9.772 

10.280 

9.958 

10.330 

9.938 

10.450 

10.012 

10.634 

10.124 

10.482 

9.870 

10.358 

9.796 

10.390 

9.804 

10.240 

9.586 

9.990 

9.364 

9.628 

8.990 

9.426 

8.806 

9.098 

8.492 

8.948 

8.244 

8.756 

8.098 

8.834 

8.202 

8.746 

7.976 

8.688 

7.968 

8.154 

7.544 

8.070 

7.444 

7.766 

7.188 

7.174 

6.684 

7.186 

6.676 

6.994 

6.388 

6.642 

6.132 

6.528 

6.036 

6.350 

5.926 

6.214 

5.750 

5.864 

5.502 

5.664 

5.224 

5.286 

4.940 

5.340 

4.960 
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TABLE  B.4:  Four- Variable  Four- Valued  Average  Output  Product  Terms 


Number  of 
Input  Terms 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood 

Decoupling 

3 

3.136 

3.018 

3.002 

4 

4.252 

4.026 

3.974 

5 

5.248 

5.028 

4.978 

6 

6.546 

6.070 

6.018 

7 

7.690 

7.168 

7.054 

8 

8.786 

8.178 

8.066 

9 

9.988 

9.172 

8.970 

10 

11.172 

10.246 

10.014 

11 

12.354 

11.330 

10.960 

12 

13.600 

12.374 

12.074 

13 

14.676 

13.502 

13.066 

14 

16.096 

14.512 

14.144 

15 

16.954 

15.366 

15.056 

16 

18.028 

16.348 

16.230 

17 

19.144 

17.570 

17.036 

18 

20.050 

18.648 

17.866 

19 

21.072 

19.460 

18.862 

20 

21.930 

20.590 

19.850 

21 

23.074 

21.396 

20.576 

22 

23.682 

22.258 

21.456 

23 

24.456 

22.910 

21.954 

24 

25.258 

23.590 

22.770 

25 

25.528 

24.552 

23.220 

26 

25.740 

24.788 

23.654 

27 

26.656 

25.412 

24.374 

28 

26.940 

26.190 

24.950 

29 

27.786 

26.668 

25.410 

30 

27.672 

26.612 

25.454 

31 

28.032 

27.246 

25.854 

32 

28.458 

27.550 

26.318 

33 

28.368 

27.840 

26.454 

34 

28.768 

27.934 

26.344 

35 

28.752 

28.170 

26.498 

48 


APPENDIX  C:  AVERAGE  COMPUTATION 
TIME  FOR  MVL  MINIMIZATION 


TABLE  C.l:  Two- Variable  Four- Valued  Average  Computation  Time 


Number  of 
Input  Terms 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood 

Decoupling 

3 

0.03309 

0.06156 

0.06336 

4 

0.05517 

0.10724 

0.11221 

5 

0.08384 

0.16931 

0.17947 

6 

0.11167 

0.23128 

0.24318 

7 

0.14140 

0.29612 

0.31576 

8 

0.K  ‘0 

0.37093 

0.38907 

9 

0.19713 

0.43512 

0.46649 

10 

0.22143 

0.48732 

0.52141 

11 

0.24440 

0.55478 

0.59490 

12 

0.27069 

0.61506 

0.65662 

13 

0.28888 

0.66258 

0.71877 

14 

0.31570 

0.71373 

0.77805 

15 

0.33116 

0.76289 

0.82645 

16 

0.34269 

0.79065 

0.85818 

49 


TABLE  C.2:  Two- Variable  Five- Valued  Average  Computation  Time 


Number  of 
Input  Terms 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood 

Decoupling 

3 

0.05604 

0.10954 

0.11110 

4 

0.09554 

0.20432 

0.20372 

5 

0.15442 

0.33730 

0.33034 

6 

0.22760 

0.50900 

0.50168 

7 

0.29088 

0.67988 

0.66278 

8 

0.36180 

0.86232 

0.84642 

9 

0.42612 

1.08536 

1.05530 

10 

0.50258 

1.24780 

1.23926 

11 

0.57016 

1.44212 

1.41004 

12 

0.62940 

1.64508 

1.60892 

13 

0.70374 

1.83080 

1.79734 

14 

0.76936 

2.07240 

1.98956 

15 

0.84158 

2.23954 

2.14072 

16 

0.88608 

2.36738 

2.27968 

17 

0.93786 

2.50298 

2.45116 

18 

0.99650 

2.63846 

2.62578 

19 

1.03754 

2.72054 

2.74588 

20 

1.07786 

2.84114 

2.86624 

21 

1.13426 

2.94840 

2.98710 

22 

1.15298 

2.99784 

3.08048 

23 

1.18984 

3.04716 

3.10488 

24 

1.23294 

3.10796 

3.25222 

25 

1.24788 

3.24612 

3.37128 

50 


Number  of 
Input  Terms 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


and  Armstrong 

Dueck  and  Miller 

Neighborhocx 

Decoupling 

0.11096 

0.25190 

0.23496 

0.17684 

0.41760 

0.39676 

0.29150 

0.67308 

0.66464 

0.43990 

1.04038 

1.02516 

0.65294 

1.49640 

1.47802 

0.86260 

2.09008 

2.07168 

1.07620 

2.61164 

2.59776 

1.36058 

3.29064 

3.24702 

1.63312 

3.97860 

3.93274 

1.97194 

4.71946 

4.73014 

2.24472 

5.53536 

5.50418 

2.61166 

6.50490 

6.37408 

3.00044 

7.56094 

7.31712 

3.40472 

8.50130 

8.33298 

3.80686 

9.56092 

9.21312 

4.22300 

10.6681 

10.1814 

4.62690 

11.9961 

11.2720 

5.10534 

12.9534 

12.0204 

5.51390 

14.0393 

3.1130 

5.91520 

15.2872 

1.0559 

6.35928 

16.1587 

5.1375 

6.69882 

17.5192 

15.9774 

7.17424 

18.2629 

17.0413 

7.60032 

19.5004 

17.6663 

7.77850 

20.4014 

18.4309 

8.32284 

21.5390 

19.5850 

8.95526 

22.5366 

20.7064 

9.04034 

23.3790 

21.0932 

9.11648 

24.1279 

21.9608 

9.48198 

24.2701 

22.3971 

9.88210 

25.3621 

23.1439 

10.1734 

26.4659 

24.2141 

10.3463 

26.5370 

24.3088 

10.8508 

27.1308 

25.6463 

11.0079 

28.1984 

25.9083 

11.2138 

28.6137 

25.8980 

11.5521 

29.2124 

27.1013 

11.7711 

30.3718 

28.0957 

11.8719 

30.1055 

28.1734 

12.4192 

30.7879 

29.1146 

12.5431 

31.1774 

29.4338 

12.3511 

30.4067 

29.0054 

12.6780  51 

31.7510 

30.1296 

TABLE  C.4:  Four- Variable  Four- Valued  Average  Computation  Time 


Number  of 
Input  Terms 

Pomper  and  Armstrong 

Dueck  and  Miller 

Neighborhood 

Decoupling 

3 

0.35778 

1.01774 

0.79966 

4 

0.57068 

1.64304 

1.33264 

5 

0.89762 

2.74510 

2.17928 

6 

7 

1.26432 

1.69976 

3.97440 

5.14632 

3.24550 

4.58174 

8 

2.52112 

7.31790 

6.54000 

9 

3.25296 

9.42768 

8.51408 

10 

4.23814 

12.6106 

11.3127 

11 

5.21906 

15.5529 

13.8771 

12 

6.31822 

18.5009 

17.0987 

13 

7.86472 

22.1662 

20.9875 

14 

9.35484 

25.1544 

24.6584 

15 

10.9688 

29.2464 

29.0192 

16 

13.0702 

34.6918 

34.6770 

17 

15.9150 

39.5156 

40.0190 

18 

17.5832 

45.6882 

46.5026 

19 

20.8329 

52.2132 

54.0485 

20 

22.9053 

58.4218 

59.8628 

21 

26.1380 

65.7590 

67.4621 

22 

30.5844 

74.2779 

76.7863 

23 

34.5081 

82.4673 

86.0243 

24 

39.0754 

92.3658 

94.9414 

25 

42.8344 

104.355 

100.501 

26 

47.3552 

111.323 

106.223 

27 

53.0040 

118.979 

113.879 

28 

57.5171 

125.230 

119.541 

29 

60.7629 

136.164 

129.975 

30 

63.8871 

146.791 

140.376 

31 

70.0774 

160.949 

154.421 

32 

73.9980 

172.364 

165.895 

33 

78.4594 

185.003 

177.013 

34 

84.8676 

200.174 

191.966 

35 

92.5689 

215.201 

207.516 
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